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In this paper we present a paradigmatic tight-binding model for single-layer as well as for multi- 
layered semiconducting M0S2 and similar transition metal dichalcogenides. We show that the elec- 
tronic properties of multilayer systems can be reproduced in terms of a tight-binding modelling 
of the single-layer hopping terms by simply adding the proper interlayer hoppings ruled by the 
chalcogenide atoms. We show that such tight-binding model permits to understand and control 
in a natural way the transition between a direct-gap band structure, in single-layer systems, to an 
indirect gap in multilayer compounds in terms of a momentum/orbital selective interlayer splitting 
of the relevant valence and conduction bands. The model represents also a suitable playground to 
investigate in an analytical way strain and finite-size effects. 
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I. INTRODUCTION 

The isolation of flakes of single-layer and few-layer 
graphenei"— has triggered a huge burst of interest on two- 
dimensional layered materials because of their structural 
and electronic properties. Due to its huge electronic mo- 
bility, graphene has been in the last years the main focus 
of the research in the field. However, a drawback in en- 
gineering graphene-based electronic device is the absence 
of a gap in the monolayer samples, and the difficulty in 
opening a gap in multilayer systems without affecting 
the mobility. As an alternative route, recent research is 
exploring the idea of multilayered heterostructures built 
up from interfacing different twodimensional materials.— 
Along this perspective, semiconducting dichalcogenides 
such as M0S2, MoSe2, WS2, etc. are promising com- 
pounds since they can be easily exfoliated and present a 
suitable small gap both in single-layer and in few-layer 
samples. Quite interestingly, in few-layer M0S2 the size 
and the nature of the gap depends on the number N 
of M0S2 layers, with a transition between a direct gap 
in monolayer [N = 1) compounds to a smaller indirect 
gap for N > 2^^— In addition, the electronic properties 
appear to be highly sensitive to the external pressure 
and strain, which affect the insulating gap and, under 
particular conditions, can also induce a insulator/metal 
transition)^"— Another intriguing feature of these mate- 
rials is the strong entanglement between the spin and the 
orbital/valley degrees of freedom, which permits, for in- 
stance, to manipulate spins by means of circularly polar- 
ized lighti^"— Moreover, in MoSe2, a transition between 
a direct to an indirect gap was observed as a function of 
temperaturei^ 

On the theoretical level, the description of its low- 
energy electronic properties is enormously facilitated by 
the the availability of a paradigmatic Hamiltonian model 
for the single-layer in terms of few tight-binding (TB) 
parameters^Si^ (actually only one, the nearest neigh- 
bors carbon-carbon hopping 70, in the simplest case)i^ 



The well-known Dirac equation can thus be derived from 
that as a low-energy expansion. Crucial to the develop- 
ment of the theoretical analysis in graphene is also the 
fact that model Hamiltonians for multilayer graphenes 
can be built using the single-layer TB description as a 
fundamental block and just adding additional interlayer 
hopping terms 1^"— Different stacking orders can be also 
easily investigated. The advantage of such tight-binding 
description with respect to first-principles calculations is 
that it provides a simple starting point for the further in- 
clusion of many-body electron-electron effects by means 
of Quantum Field Theory (QFT) techniques, as well as of 
the dynamical effects of the electron-lattice interaction. 
Tight-binding approaches can be also more convenient 
than first-principles methods such as Density Functional 
Theory (DFT) for investigating systems involving a very 
large number of atoms. Although DFT methods are cur- 
rently able to handle systems with hundreds or even thou- 
sands of atoms^i^, and have been thoroughly applied to 
large scale graphene-related problems^"—, they are still 
computationally challenging and demanding. Therefore, 
TB has been the method of choice for the study of disor- 
dered and inhomogeneous systems^"— materials nanos- 
tructured in large scales (nanoribbons, ripples)^"— or in 
twisted multilayer materialsi^"— 

While much of the theoretical work of graphenic ma- 
terials has been based on tight-binding-like approaches, 
the electronic properties of single-layer and few-layer 
dichalcogenides have been so far mainly investigated by 
means of DFT calculations <^i^i2r2SSi2^, despite early 
work in non-orthogonal tight binding models for transi- 
tion metal dichalcogenides 1^ Few simplified low-energy 
Hamiltonian models has been presented for these ma- 
terials, whose validity is however restricted to the spe- 
cific case of single-layer systems. An effective low-energy 
model was for instance introduced in Refs. |24||97| to dis- 
cuss the spin/orbital/valley coupling at the K point. Be- 
ing limited to the vicinity of the K point, such model can- 
not be easily generalized to the multilayer case where the 



gap is indirect with valence and conduction edges located 
far from the K point. An effective lattice TB Hamilto- 
nian was on the other hand proposed in Refs. [9a, valid 
in principle in the whole Brillouin zone. However, the 
band structure of the single-layer lacks the characteristic 
second niininium in the conduction band (see later dis- 
cussion) that will become the effective conduction edge 
in multilayer systems, so that also in this case the gen- 
eralization to the multilayer compounds is doubtful. In 
addition, the use of an overlap matrix makes the pro- 
posed Hamiltonian unsuitable for a straightforward use 
as a basis for QFT analyses. This is also the case for a re- 
cent model proposed in Ref. |99|, where the large number 
(ninetysix) of free fitting parameters and the presence of 
overlap matrix make such model inappropriate for prac- 
tical use within the context of Quantum Field Theory. 

In this paper we present a suitable tight-binding model 
for the dichalcogenides valid both in the single-layer 
case and in the multilayer one. Using a Slater-Koster 
approachjiSS and focusing on M0S2 as a representative 
case, we analyze the orbital character of the electronic 
states at the relevant high-symmetry points. Within this 
context wc show that the transition from a direct gap 
to an indirect gap in M0S2 as a function of the number 
of layers can be understood and reproduced in a natural 
way as a consequence of a momentum/orbital selective 
interlayer splitting of the main relevant energy levels. In 
particular, we show that the Pz orbital of the S atoms 
plays a pivotal role in such transition and it cannot be 
neglected in reliable tight-binding models aimed to de- 
scribe single-layer as well as multi-layer systems. The 
tight-binding description here introduced can represent 
thus the paradigmatic model for the analysis of the elec- 
tronic properties in multilayer systems in terms of intra- 
layer ligands plus a finite number of interlayer hopping 
terms. Such tight-binding model, within the context of 
the Slater-Koster approach, provides also a suitable tool 
to include in an analytical and intuitive way effects of 
pressure/strain by means of the modulation of the in- 
teratomic distances. The present analysis defines, in ad- 
dition, the minimum constraints that the model has to 
fulfill to guarantee a correct description of the band struc- 
ture of multi-layer compounds. 

The paper is structured as follows: in Section |lT] we 
present DFT calculations for single-layer and multi-layer 
(bulk) M0S2, which will be here used as a reference for 
the construction of a tight-binding model. In Section 
mil we describe the minimum tight-binding model for the 
single-layer case needed to reproduce the fundamental 
electronic properties and the necessary orbital content. 
The decomposition of the Hamiltonian in blocks and the 
specific orbital character at the high-symmetry points is 
discussed. The extension of the tight-binding model to 
the bulk case, taken as representative of multilayer com- 
pounds, is addressed in Section HVl We pay special at- 
tention to reveal the microscopic origin of the change 
between a direct-gap to indirect-gap band structure. In 
Section |V] we summarize the implications of our analysis 
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FIG. 1: (a) Model of the atomic structure of M0S2. The bulk 
compound has a 2H-M0S2 structure with two M0S2 layers 
per unit cell, each layer being built up from a trigonal prism 
coordination unit. The small green rectangle represents the 
unit cell of a monolayer of M0S2, which is doubled (red ex- 
tension) in the bulk crystal, (b) Detail of the trigonal prisms 
for the two layers in the bulk compound, showing the lattice 
constants and the definition of the structural angles used in 
the text. 



in the building of a reliable tight-binding model, and we 
provide a possible set of tight-binding parameters for the 
single-layer and multilayer case. 



II. DFT CALCULATIONS AND ORBITAL 
CHARACTER 

In the construction of a reliable TB model for semi- 
conducting dichalcogenides we will be guided by first- 
principles DFT calculations taht will provide the refer- 
ence on which to calibrate the TB model. We will fo- 
cus here on M0S2 as a representative case, although we 
have performed first-principle calculations for compari- 
son also on WS2. The differences in the electronic struc- 
ture and in the orbital character of these two compounds 
are, however, minimal and they do not involve any differ- 
ent physics. The structure of single-layer and multilayer 
M0S2 is depicted in Fig. [T] 

The basic unit block is composed of an inner layer of 
Mo atoms on a triangular lattice sandwiched between 
two layers of S atoms lying on the triangular net of al- 
ternating hollow sites. Following standard notations^ 
we denote a as the distance between nearest neighbor in- 
planc Mo-Mo and S-S distances, h as the nearest neighbor 
Mo-S distance and u as the distance between the Mo and 
S planes. The M0S2 crystal forms an almost perfect trig- 
onal prism structure with 6 and u very close to the their 
ideal values b ~ \Jl jYla and u ~ a/2. In our DFT cal- 
culations, we use experimental values for bulk M0S2,— 
namely a = 3.16 A, u = 1.586 A, and, in bulk systems, 
a distance between Mo planes as c' = 6.14 A, with a 
lattice constant in the 2H-M0S2 structure of c = 2c'. 




FIG. 2: Band structure and orbital character of single-layer 
M0S2. The top left panel shows the full band structure while, 
in the other panels, the thickness of the bands represents the 
orbital weight, where the d-character {d2 — d^2_y2, dxy, di = 
dxz, dyz, do — d^-,2_^2) refers to the Mo atom 4d orbitals, 



while the p-character {p^y 
sulfur. 



Px, Py) refers to 2p orbitals of 



The in-plane Brillouin zone is thus characterized by the 
high-symmetry points F = (0,0), K= 47r/3a(l,0), and 
M= 47r/3a(0, •\/3/2). DFT calculations are done using 
the Siesta codci ^^i"^^ We use the exchange-correlation 
potential of Ceperly-AlderiiSi as parametrized by Perdew 
and ZungerjiS^ We use also a split-valence doublc-C basis 
set including polarization functionsJ^ The energy cutoff 
and the Brillouin zone sampling were chosen to converge 
the total energy. 

The electronic dispersion for the single-layer M0S2 is 
nowadays well known. We will only focus on the block of 
bands containing the first four conduction bands and by 
the first seven valence bands, in an energy window of from 
-7 to 5 eV around the Fermi level. Our DFT calculations 
are shown in Fig. [2J where we show the orbital character 
of each band. We use here the shorthand notation ^2 
to denote Mo 4(1^2 _y2, 4 da; y orbitals; di for the Mo Adxz, 
Adyz orbitals; do for the Mo 4d3^2_r2 orbital; p^y (or 
simply p) to denote the S 3pa;, ipy orbitals; and pz (or 
simply z) for the S Sp^orbital. The four conduction bands 
and the seven valence bands are mainly constituted by 
the five 4d orbitals of Mo and the six (three for each 
layer) 3p orbitals of S, which sum up to the 93 % of the 
total orbital weight of these bands. 

A special role in the electronic properties of these ma- 
terials is played by the electronic states labeled as (A)- 
(D) and marked with black bullets in Fig. [5J A detailed 
analysis of the orbital character of each energy level at 
the main high-symmetry points of the Brillouin zone, as 
calculated by DFT, is provided in Table IH We can notice 
that an accurate description of the conduction and va- 
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* Double-degenerate level 

TABLE I: Energy levels and orbital content of single-layer 
M0S3 evaluated by DFT calculations. We report here the 
first two main orbital characters belonging to the blocks Mo- 
4d and S-3p, while the following column shows the remaining 
character not belonging to these orbital group. Also show is 
the association of each level with the corresponding eigenvalue 
of the tight-binding model and the symmetry with respect to 
the z -^ —z inversion (E=even, 0=odd). The label Eai3,± 
in the last column denotes the orbital character of the TB 
eigenstate, with a, P — p, 2:,da, di, do, where p = Px,Py, z = 
Pz, da = dx2_y2,dxy, di = dxz, dyz, do = d3^2_^2. The index 
± denotes the higher energy [(-I-) = antibonding] and the 
lower energy [( — ) = bonding]. 



lence band edges (A)-(B) at the K point involves at least 
the Mo orbitals d^z^-r^, da;2_j,2, dxy, and the S orbitals 
Px, Py Along this perspective, a 5-band tight-binding 
model, restricted to the subset of these orbitals, was pre- 
sented in Ref. [9^, whereas even the S 3p orbitals were 
furthermore omitted in Ref. \24 . 

The failure of this latter orbital restriction for a more 
comprehensive description is however pointed out when 
analyzing other relevant high-symmetry Brillouin points. 
In particular, concerning the valence band, we can notice 
a second maximum at the F point, labeled as (C) in Fig. 
[51 just 42 meV below the real band edge at the K point 
and with main do-pz orbital character. The relevance of 
this secondary band extreme is evident in the multilayer 
compounds {N > 2), where such maximum at F increases 
its energy to become the effective band edgci^i^ 

The band structure with the orbital character for the 
bulk [N = 00) case, representative of the multilayer case, 
is shown in Fig. [S] A similar change of the topol- 
ogy of the band edge occurs in the conduction band. 
Here a secondary minimum, labeled as (D) in Fig. [^ at 
Q = 47r/3a(l/2, 0), midway along the F-K cut, is present 
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FIG. 3: Band structure and orbital character for bulk 2H- 
M0S2. Labels similar as in Fig. [S] 



in the single-layer compounds. Such minimum however 
moves down in energy in multilayer systems to become 



the effective conduction band edgci^ Even in this case, 
a relevant pz component is involved in the orbital char- 
acter of this electronic state. The topological changes of 
the location of the band edges in the Brillouin zone are 
responsible for the observed switch from a direct to an 
indirect gap in multilayer samples. As we will see, thus, 
the inclusion of the Pz orbitals in the full tight-binding 
Hamiltonian is not only desirable for a more complete 
description, but it is also unavoidable to understand the 
evolution of the band structure as a function of the num- 
ber of layers. 



III. TIGHT-BINDING DESCRIPTION OF THE 
SINGLE-LAYER 

The aim of this section is to define a tight-binding 
model for the single-layer which will be straightforwardly 
generalizable to the multilayer case by adding the appro- 
priate interlayer hopping. We will show that, to this 
purpose, all the Ad Mo orbitals and the 3p S orbitals 
are needed to be taken into account. Considering that 
the unit cell contains two S atoms, we define the Hilbert 
space by means of the 11-fold vector: 



J 



(pI, 



'i,x,t^Pi,y,vPi,z,t^ "'i,3z^-r^^ 



' ^i,xz^ VJ/z 



^Pi,x,b^Pi,y,b'Pi,z,b)^ 



(1) 



r 



where di^a creates an electron in the orbital a of the Mo 
atom in the i-unit cell, Pi^cf,t creates an electron in the 
orbital a of the top (<) layer atom S in the i-unit cell, and 
Pi,a,b creates an electron in the orbital a of the bottom 
(&) layer atom S in the i-unit cell. 



Once the Hilbert space has been introduced, the tight- 
binding model is defined by the hopping integrals be- 
tween the different orbitals, described, in the framework 
of a Slater-Koster description, in terms of a, tt and 5 
ligandsJ^ In order to provide a tight-binding model as 
a suitable basis for the inclusion of many-body effects by 
means of diagrammatic techniques, we assume that the 
basis orbitals are orthonormal, so that the overlap matrix 
is the unit matrix. A preliminary analysis based on the 
interatomic distance can be useful to identify the most 
relevant hopping processes. In particular, these are ex- 
pected to be the ones between nearest neighbor Mo-S (in- 
teratomic distances b = 2.41 A) and between the nearest 
neighbor in-plane Mo-Mo and between the nearest neigh- 
bor in-plane and out-of-plane S-S atoms (interatomic dis- 
tance a = 3.16 A). Further distant atomic bonds, in 
single-layer systems, start from hopping between second 
nearest neighbor Mo-S atoms, with interatomic distance 
3.98 A, and they will be here discarded. 



All the hopping processes of the relevant pair of neigh- 
bors are described in terms of the Slater-Koster param- 
eters, respectively Vpda, Vpd-^ (Mo-S bonds), Vdda, Vddn, 
VddS (Mo- Mo bonds), and Vppc, Vpp-K (S-S bonds). Addi- 
tional relevant parameters are the crystal fields Ag, Ai, 
A2, Ap, Az, describing respectively the atomic level the 
I = ((i322_r2), the I = 1 {dxz, dyz), the I = 2 {dx'^-y^, 
dxy) Mo orbitals, the in- plane {px, Py) S orbitals and of 
the out-of-plane pz S orbitals. We end up with a total of 
12 tight-binding parameters to be determined, namely: 

Aq, Ai, A2, Ap, A2, Vdda, VddiT, VddS, Vppa, Vpp.^, Vpda, 

In the orbital basis of Eq. ([T]), we can write thus the 
tight-binding Hamiltonian in the form: 



H 



E' 



>k-ffk0k, 



(2) 



where (/)k is the Fourier transform of 4>i in momentum 
space. The Hamiltonian matrix can be written (we drop 
for simplicity from now on the index k) as: 



H 



J^pt,pt ^d,pt ^pt,pb 

Hd,pt Hd.d Hd,pb 

^pb,pb ^d,pb ^pb,pb 



(3) 



where H, 



pb,ph 



Hpt^pt describes the in-plane hopping in 



the top and bottom S layer, namely, 



H, 



pb,pb 
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Hd^d the in-plane hopping in the middle Mo layer, 
namely, 
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(5,) 



xz/yz 



tiyz/yz I 



Hpt.pb the vertical hopping between S orbitals in the top 
and bottom layer. 



H, 



pt,pb 
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K, 



(6) 



and Hd.pt, Hd^pb the hopping between Mo and S atoms 
in the top and bottom planes, respectively: 
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yz/x 



Here and in the following, for the sake of compact- 
ness, we use the shorthand notation 3z^ — r^ ^ z^ and 
x^ — y^ ^ x'^. An explicit expression for the differ- 
ent Hamiltonian matrix elements in terms of the Slater- 
Koster tight-binding parameters can be provided follow- 
ing the seminal work by Doran et al. (Ref. Il04f ) and it 
is reported for completeness in Appendix lAl 

Eqs. (HI)-® define our tight-binding model in terms 
of a 11 X 11 Hamiltonian H which can be now explicitly 
solved to get eigenvalues and eigenvectors in the whole 
Brillouin zone or along the main axes of high symmetry. 
It is now an appealing task to associate each DFT energy 
level with the Hamiltonian eigenvalues, whose eigenvec- 
tors will shed light on the properties of the electronic 
states. Along this line, we are facilitated by symmetry 
arguments which permit, in the monolayer compounds, 
to decoupled the 11 x 11 Hamiltonian in Eq. ([3]), in two 
main blocks, with different symmetry with respect to the 
mirror inversion z — > — z.ii24 This task is accomplished by 
introducing a symmetric and antisymmetric linear combi- 
nation of the p orbital of the S atoms on the top/bottom 
layers. More explicitly, we use the basis vector 



(4 



fc,3z2 



-<: 



y2^"-k,xy^Pk,x,S'Pk,y,S^Pk,z,A^"- 



k,xz 



^"'k,yz^Pk,x,A^Pk,y,A^Pk,z,Sh 



(9) 



r 



where pl^^^g = [pl^^^ + pi,„, J/V2, pl^^^ = bl,a,t 



pI 



^)/-\/2- Note that our basis di ffers slightly with re- 
spect to the one employed in Ref. Il04l because we have 
introduced explicitly the proper normalization factors to 
make it unitary. In this basis we can write thus 



high-symmetry points of the Brillouin zone. Most im- 
portant are the K and the F points, which are strictly 
associated with the direct and indirect gap in monolayer 
and multilayered compounds. 



H = 



He 
Ho 



(10) 



where H^, is a 6 x 6 block with even (E) symmetry with 
respect to the mirror inversion 2; — > — z, and Hq a 5 x 5 
block with odd (O) symmetry. We should remark how- 
ever that such decoupling holds true only in the single- 
layer case and only in the absence of a z-axis electric 
field, as it can be induced by substrates or under gating 
conditions. In the construction of a tight-binding model 
that could permit a direct generalization to the multilayer 
case, the interaction between the band blocks with even 
and odd symmetry should be thus explicitly retained. 

The association between DFT energy levels and tight- 
binding eigenstates is now further simplified on specific 



A. r point 



We present here a detailed analysis of the eigenstates 
and their orbital character at the F point. For the sake of 
simplicity, we discuss separately the blocks with even and 
odd symmetry with respect to the inversion z ^ —z. The 
identification of the DFT levels with the tight-binding 
eigenstates is facilitated by the possibility of decomposing 
the full Hamiltonian in smaller blocks, with typical size 
2x2 (dimers) or 1 x 1 (monomers). In particular, the 
6x6 block with even symmetry can be decomposed (see 



Appendix [B] for details) as: 



He{T) 





-ffpd.(r) 
^pd.(r) 



(11) 



Here each matrix, -f/pdj, H^da represents a 2 x 2 block 
where the indices describe the orbital character of the 
dimer. In particular, Hp^^ involves only d2 = d^2_y2,dxy 
Mo-orbitals and Px,Py S-orbitals, whereas Hzdo involves 
only the do = d^^^^r^ Mo-orbital and the pz S-orbital. 
As it is evident in ((TTj) . the block Hpd2 appears twice 
and it is thus double degenerate. Similarly, we have 



^o(r) 



Wi(r 
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HpdA^) 
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(12) 



where the doubly degenerate block Hpdi involves only 
di = dxz, dyz Mo-orbitals and Px,Py S-orbitals, while F^ 
is a 1 X 1 block (monomer) with pure character pz . 

It is also interesting to give a closer look at the inner 
structure of a generic Hamiltonian sub-block. Consider- 
ing for instance Hzdo ^^ ^^ example, we can write 



H^doir) 



To v2Tzdo 

V2rzdo rE 



(13) 



where Fg is an energy level with pure Mo do orbital char- 
acter and Ff an energy level with pure S Pz orbital char- 
acter. The off-diagonal term ^/2Tzdo acts thus here as a 
"hybridization" , mixing the pure orbital character of Fg 
and Ff . The suffix "E" here reminds that the level F? 
belongs to the even symmetry block, and it is useful to 
distinguish this state from a similar one with odd symme- 
try (and different energy). Keeping Hzda as an example, 
the eigenvalues of a generic 2x2 block can be obtained 



analytically: 



E.do,±0^) = 



± 



To-Tf 



-2rL,(i4) 



The explicit expressions of Fq and Ta/3 in terms of the 
Slater-Koster tight-binding parameters is reported in Ap- 
pendix |X1 

It is interesting to note that the diagonal terms Fq, 
(a = do,di,d2,p, z) are purely determined by the crys- 
tal fields Aq, and by the tight-binding parameters Vdda, 
Vddn, VddS, Vppa, VppT,, connecting Mo- Mo and S-S atoms, 
whereas the hybridization off-diagonal terms V^p depend 
exclusively on the Mo-S nearest neighbor hopping Vpda, 

ypd-K- 

A careful comparison between the orbital character of 
each eigenvector with the DFT results permits now to 
identify in an unambiguous way each DFT energy level 
with its analytical tight-binding counterpart. Such asso- 
ciation is reported in Table U where also the even/odd 
symmetry inversion is considered. 

The use of the present analysis to characterize the 
properties of the multilayer M0S2 will be discussed in 
Section HVl 



B. K point 

A crucial role in the properties of semiconducting 
dichalcogenides is played by the K point in the Brillouin 
zone, where the direct semiconducting gap occurs in the 
single-layer systems. The detailed analysis of the elec- 
tronic spectrum is also favored here by the possibility 
of reducing the full 11x11 Hamiltonian in smaller sub- 
blocks. This feature is, however, less evident than at the 
F point. The even and odd components of the Hamilto- 
nian take the form: 



He{K) 



Ho{K) 



Ko -i^/2Kpao V^Kpd, 




K2 i^Kpd, V2Kpd, V2Kzd2 


K2 -V2Kpd, iV2Kpd2 -iV2Kzd 


iV2Kpdo -iV2Kpd, -V2Kpd, K^ 


V2Kpa„ V2Kpd, -iV2Kpd, K^ 


V2Kzd2 iV2Kzd2 Kf 


/ ^1 V2Kpd, ~iV2Kpd, -iV2Kzd, \ 






Ki -iV2Kpd, ~V2Kpd, V2Kzd, 




= 


V2Kpd, iV2Kpd, A'O 






iV2Kpd, -V2Kpd, K^ 






\iV2Kzd, V2Kzd, ATf y 





(15) 



(16) 



As for the F point, also here the upper labels (/i =E, 
O) in K^ {^ =E, O) express the symmetry of the state 



corresponding to the energy level K^ with respect to 



the z — >■ — z inversion. The electronic properties of the 
Hamihonian at the K point look more transparent by 

I 



introducing a different "chiral" base: 
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(^fc,3z2 _r2 ' "fc,L2 ' "'kM2 ' Pk,L,S^ Pk,R,S ^ Pk, 



r 



,A'"fc.Ll'"fc. 



Rl 



•Pk,L,A^Pk,R,A^Pk,z,s)^ 



(17) 



where (ifc,L2 = {dk^x-2_y-2-idk^xy) / \/^, rffc,_R2 == {dk,x-^~y-^ + 
idk,xy)/V^, dk,Li = idk,xz~idk,yz)/V2, dk,Ri = idk,xz + 

idk,yz)/V^, Pk,L,S = iPk,x,S - iPk,y,s)/V^, Pk,R,S = 
{Pk,x,S + iPk^y^s)/"/^, Pk,L,A = ij>k,x,A - ipk,y,A)/\/^, 
Pk,R,A = {Pk,x,A + iPk,y,A)l^- 

In this basis, the Hamiltonian matrix can be also di- 
vided in smaller sub-blocks (see Appendix [BJ as: 



ii^{K) 



and 



-ffo = 



^pdo(i^) 

H,i^{K) 

%ciAi<) 



HpdAK) 

H.aAK) 

K^ 



(18) 



(19) 



As it is evident from the labels, each sub-block is also 
here a 2 x 2 dimer, apart from the term K^ which is a 
1x1 block (monomer) with pure Px^Py character. The 
association between the DFT energy levels and the tight- 
binding eigenstates is reported also for the K point in 
Table i 



C. Q point 

As discussed above, another special point determining 
the electronic properties of M0S2 is the Q point, halfway 
between the T and K points in the Brillouin zone, where 
the conduction band, in the single-layer system, has a 
secondary minimum in addition to the absolute one at 
the K point. Unfortunately, not being a point of high- 
symmetry, the tight-binding Hamiltonian cannot be de- 
composed in this case in simpler smaller blocks. Each 
energy eigenvalue will contain thus a finite component of 
all the Mo and S orbitals. In particular, focusing on the 
secondary minimum in Q, DFT calculations give 46 % 
^2, 24 % Px/y, 11% Pz and 9 % do- The orbital content 
of this level will play a crucial role in determining the 
band structure of multilayer compounds. 



D. Orbital constraints for a tight-binding model 

After having investigated in detail the orbital contents 
of each eigenstate at the high-symmetry points, and hav- 
ing identified them with the corresponding DFT energy 



levels, we can now employ such analysis to assess the 
basilar conditions that a tight-binding model must fulfill 
and to elucidate the physical consequences. 

A first interesting issue is about the minimum number 
of orbitals needed to be taken into account in a tight- 
binding model for a robust description of the electronic 
properties of these materials. A proper answer to such 
issue is, of course, different if referred to single-layer or 
multilayer compounds. For the moment we will focus 
only on the single-layer case but we will underline on the 
way the relevant features that will be needed to take into 
account in multi-layer systems. 

In single-layer case, focusing only on the band edges 
determined by the states (A) and (B) at the K point, 
we can identify them with the eigenstates Epdo.+{K), 
Epd2.+ {K), respectively, with a dominant Mo 'id char- 
acter and a marginal S Px/y component, as we show be- 
low. It is thus tempting to define a reduced 3-band tight- 
binding model, keeping only the Mo 4d3j,2_^2, 4d^2_j^2, 
Adxy orbitals with dominant character and disregarding 
the S Px^ Py orbitals, with a small marginal weight. A 
similar phenomenological model was proposed in Ref. |24| . 
However, the full microscopic description here exposed 
permits to point out the inconsistency of such a model. 
This can be shown by looking at Eq. (fTS)) . The band gap 
at K in the full tight- binding model including S Px, Py 
orbitals is determined by the upper eigenstate of ffpdo 1 



Epda,+ {K) - 



Ko + K^ 




Ko - Kf 



4i^|.o' (20) 



and the upper eigenstate of Hpd2 , 

Epd2.+ {K) 



K2 + K^ 



'K2-K^ 



" ^ +SKl,^, (21) 



both with main Mo 4d character, while the eigenstate 
K2 + K^ 



Ezd^AK) 




K2-Kf 



4i^L' (22) 





K 



K 



FIG. 4: Schematic band structure close to the K point for the 
valence and conduction bands: (a) including S Px, Py orbitals; 
(b) omitting S Px, Py orbitals. 



also with dominant Mo 4d character, but belonging to 
the block Hzd2i lies at higher energy (see table |I]). The 
3-band model retaining only the do, d2 orbitals is equiva- 
lent to switch off the hybridization terms Kpdg , Kpd2 j 
K^d2, ruled by Vpda, Vpd-^, so that Epdo.+{K) = Kq, 
Epd2,+{K) = Ezd2,+ {K) = K2- In this context the level 
Ezd2,+ {K) becomes degenerate with Epd2,+ {K). This de- 
generacy is not accidental but it reflects the fact that the 
elementary excitations of the ^2 states, in this simplified 
model, are described by a Dirac spectrum, as sketched 
in Fig. m As a consequence, no direct gap can be possi- 
bly established in this framework. It is worth to mention 
that a spin-orbit coupling can certainly split the Dirac 
cone to produce a direct gap at the K point, but it would 
not explain in any case the direct gap observed in the 
DFT calculations without spin-orbit coupling. 

We should also mention that, in the same reduced 3- 
band model keeping only the do and d2 Mo orbitals, the 
secondary maximum (C) of the valence band would have 
a pure do orbital character. As we are going to see in the 
discussion concerning the multilayer samples, this would 
have important consequences on the construction of a 
proper tight-binding model. 

A final consideration concerns the orbital character 
of the valence band edge, Epd2,+ {K)- This state is as- 
sociated with the third 2x2 block of (fT8|) and it re- 
sults from the hybridization of the chiral state df^ji2 = 
{dk.j.2_y2 +idk.xy)/V^ of the Mo d orbitals with the chiral 
state Pk,R,s = {Pk,x,s + 'iPk,y,s)/V2 of the S p orbitals. 
The role of the chirality associated with the d orbitals, 
in the presence of a finite spin-orbit coupling, has been 
discussed in detail in relation with spin/valley selective 
probesi^"— What results from a careful tight-binding 
description is that such d-orbital chirality is indeed en- 
tangled with a corresponding chirality associated with 
the S p orbitals. The possibility of such entanglement, 
dictated by group theory, was pointed out in Ref. |27|. 

A similar feature is found for the conduction band 
edge, Epdo,+ {K). So far, this state has been assumed to 
be mainly characterized by the d^^'^-r'^, and hence with- 



out an orbital moment. However, as we can see, this is 
true only for the Mo d part, whereas the S p component 
does contain a finite chiral moment. On the other hand, 
the spin-orbit associated with the S atoms as well as with 
other chalcogenides (ex.: Se) is quite small, and taking 
into account also the small orbital S weight, the possibil- 
ity of a direct probe of such orbital moment is still to be 
explored. 



IV. BULK SYSTEM 

In the previous section we have examined in detail 
the content of the orbital character in the main high- 
symmetry points of the Brillouin zone of the single-layer 
M0S2, to provide theoretical constraints on the construc- 
tion of a suitable tight-binding model. Focusing on the 
low-energy excitations close to the direct gap at the K 
point, we have seen that a proper model must take into 
account at least the three Mo orbitals d322_^2, dx2-y2, 
dxy and the two S orbitals Px, Py On the other hand, 
our wider aim is to introduce a tight-binding model for 
the single-layer that would be the basilar ingredient for a 
tight-binding model in multilayer systems, simply adding 
the interlayer coupling. 

For the sake of simplicity we focus here on the bulk 
2H-M0S2 structure as a representative case that contains 
already all the ingredients of the physics of multilayer 
compounds. The band structure for the bulk compound 
is shown in Fig. [S] As it is known, the secondary max- 
imum (C) of valence band at the F point is shifted to 
higher energies in multilayer systems with respect to the 
single-layer case, becoming the valence band maximum. 
At the same time also the secondary minimum (D) of 
the conduction band, roughly at the Q point, is lowered 
in energy, becoming the conduction band minimum. All 
these changes result in a transition between a direct gap 
material in single-layer compounds to indirect gap sys- 
tems in the multilayer case. Although such intriguing 
feature has been discussed extensively and experimen- 
tally observed, the underlying mechanism has not been 
so far elucidated. We will show here that such topological 
transition of the band edges can be naturally explained 
within the context of a tight-binding model as a result 
of an orbital selective (and hence momentum dependent) 
band splitting induced by the interlayer hopping. 

The orbital content of the bulk band structure along 
the same high-symmetry lines as in the single-layer case 
is shown in Fig. |31 We will focus first on the K point, 
where the single-layer system has a direct gap. We note 
that the direct gap at K is hardly affected. The interlayer 
coupling produces just a very tiny splitting of the valence 
band edge Epd2,+ [K), while the conduction band edge 
Epdo,+{K) at K becomes doubly degenerate. 

Things are radically different at the F point. The anal- 
ysis of the orbital weight d^^2_^2 in Fig. [3] shows indeed 
that there is a sizable splitting of the Ezdo,+ ^) level, of 
the order of 1 eV. A bit more difficult to discern, because 



of the multi-orbital component, but still visible, is the 
splitting of the secondary minimum (D) of the conduc- 
tion band in Q. This is clearest detected by looking in 
Fig. [3] at the d2 and do characters, which belong uni- 
cally to the E block. One can thus estimate from DFT a 
splitting of this level at the Q point of ~ 1.36 eV. 

We are now going to see that all these features are 
consistent with a tight-binding construction where the 
interlayer hopping acts as an additional parameter with 
respect to the single-layer tight-binding model. From the 
tight-binding point of view, it is clear that the main pro- 
cesses to be included are the interlayer hoppings between 
the external S planes of each M0S2 block. This shows 
once more the importance of including the S p orbital 
in a reliable tight-binding model. Moreover, for geomet- 
ric reasons, one could expect that the interlayer hopping 
between the pz orbitals, pointing directly out-of-plane, 
would be dominant with respect to the interlayer hopping 
between p^, Py This qualitative argument is supported 
by the DFT results, which indeed report a big splitting 
of the Ezd„,+ 0^) level at the F point, with a 27 % of 
Pz component, but almost no splitting of the degenerate 
Epd2,+ (r) at ~ 2 eV, with 68 % component oi Px, Py 

We can quantify this situation within the tight-binding 
description by including explicitly the interlayer hopping 
between the p-orbitals of the S atoms in the outer planes 
of each M0S2 layer, with interatomic distance d = 3.49 
A (see Fig. [T]). These processes will be parametrized in 
terms of the interlayer Slater-Koster ligands Upp^, Upp.,^. 
The Hilbert space is now determined by a 22-fold vector, 
defined as: 



^I 



^l. 



~A.2l 



(23) 



where 4>\, ^ represents the basis ^ for the layer 1, and 



At 



Y-f, 2 the same quantity for the layer 2. The corresponding 
Hamiltonian, in the absence of interlayer hopping, would 
read thus: 



H, 



bulk 



Hi 

6 H2 



(24) 



where Hi , H2 refer to the intralayer Hamiltonian for the 
layer 1 and 2, respectively. 

Note that the Hamiltonian of layer 2 in the 2H-M0S2 
structure is different with respect to the one of layer 1. 
From a direct inspection we can see that the elements 
H2,a.i3{£,,il) of layer 2 are related to the corresponding 
elements of layer 1 as: 



H2.aAC,V) = Pc.PpHi^aM^,-7]), 



(25) 



where ^ = kxa/2, -q = \/3kya/2, and Pa ~ 1 if the orbital 
a has even symmetry for y ^ —y, and P„ = — 1 if it has 
odd symmetry. We note that both effects can be re- 
adsorbed in a different redefinition of the orbital basis so 
that the eigenvalues of H2 are of course the same as the 
eigenvalues of Hi . 



Taking into account the inter-layer S-S hopping terms, 
we can write thus: 



Hbu: 



Ik 



Hi Hi 

Hi H2 )' 



(26) 



where H±^ is here the interlayer hopping Hamiltonian, 
namely: 



i/i 



where C, = kzc/2 and 



/e cos C /eo sin C 
'eo 



-/j?oSinC /qCOsC 



/e = 



/o = 



'EO 



/ = 



03x3 03x3 
63x3 / 

02x2 02x3 
03X2 / 

03x2 03x3 

63x2 ii 



^x/x X I y ^x I z 
^x/y ^y/y ^y/z 
^x/z ^y/z ^z/z 



(27) 



(28) 



(29) 



(30) 



(31) 



The analytical expression of the elements Ia/i3 as func- 
tions of the Slater-Koster interlayer parameters Uppa, 
Upp-n is provided in Appendix [X] Note that, in the pres- 
ence of interlayer hopping in the bulk M0S2, we cannot 
divide anymore, for generic momentum k, the 22 x 22 
Hamiltonian in smaller blocks with even and odd sym- 
metry with respect to the change z — > — z. The analysis 
is however simplified at specific high-symmetry points of 
the Brillouin zone. In particular, for fc^ = (C = 0), we 
can easily see from (|?7|) that the block 12x12 (6x6+6x6) 
with even symmetry and the block 10 x 10 (5 x 5 + 5 x 5) 
with odd symmetry are still decoupled. 

Exploiting this feature, we can now give a closer look 
at the high-symmetry points. 



A. r point 

In Section IIIII we have seen that at the F point the 
Hamiltonian can be decomposed in 2 x 2 blocks. Par- 
ticularly important here is the block Hzda whose upper 
eigenvalue Ezdo,+ {^)i with main orbital character d^z'^-r'^ 
and a small pz component, represents the secondary max- 
imum (C) of the valence band. A first important prop- 
erty to be stressed in bulk systems is that, within this 
(Mo 4d)-|-(S 2>p) tight-binding model, the interlayer cou- 
pling at the F point does not mix any additional orbital 
character. This can be seen by noticing that the inter- 
layer matrix / is diagonal at the F point. Focusing on 



the EzdaO^) levels, we can write thus a 4 x 4 reduced 
Hamiltonian (see Appendix |B|) : 



H, 



where F, 



Uppa, UppTT between p. 





V r,. 



\ 

r,. 

To V2Tzd„ 



(32) 



represents the interlayer hopping mediated by 
orbitals belonging to the outer S 
planes on different layers. Eq. ([5^ is important because 
it shows that the qualitative idea that each energy level 
in the bulk system is just split by the interlayer hopping 
is well grounded. In particular, under the reasonable hy- 
pothesis that the interlayer hopping is much smaller than 
intralayer processes, denoting E^do.+ai^), Ezdo,+b(^) the 
two eigenvalues with primary do components, we get: 



AE. 



zdo,+ 



(r) — E;,do,+a{T) - Ezdo,+b(r) 




A similar situation is found for the other 2x2 blocks 
Hpd^iT), Hpd^iT), and the 1 x 1 block .ffz(r). Most im- 
portant, tracking the DFT levels by means of their or- 
bital content, we can note that both levels Ezdo.+ (r) and 
Ezdo,-(r) undergo a quite large splitting sa 1.2 eV, and 
the level £'z(r) a splitting sa 2.6 eV, whereas the levels 



J 
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-^p(i2(r), Hpdi{T) are almost unsplit. This observation 
strongly suggest that, as expected, the interlayer hop- 
ping between p^, py orbitals is much less effective than 
the interlayer hopping between pz . 

Similar conclusion can be drawn from the investigation 
of the energy levels at the K point, although the analysis 
is a bit more involved. 



B. K point 

The properties of the bulk system at the K point are 
dictated by the structure of the interlayer matrix / which, 
in the basis defined in Eq. ([23|) . at the K point reads: 



MK) 



iK„ 



pp 



^pz 



iK, 



pz 



— K K 

^^pp ^^pz 

Kpz 



(34) 



As discussed in detail in Appendix [B] the electronic 
structure is made more transparent by using an appro- 
priate chiral basis, which is a direct generalization of the 
one for the single-layer. We can thus write the even and 
odd parts of the resulting Hamiltonian in the form: 



Ho{K) 
where 



HpzdoAK) 

Hpzd^AK) 1(35) 

Hpd.MK) 

HpzdAK) 

HpzdAK) I (36) 

Hpd,^o{K) 



HpzdfaiK) - 



Hpd2,E{K) 



HpzdAK) 



HpduoiK) 



-2iKpd„ \ 

K^ iV2K.. 



K2 2Kzd, 



( Ko 

2iKpda 


V -iy/2Kpz 2Kzd, Kf J 

( K2 iVSKpd, \ 

K^ 2K, 

K2 



-i\/8Kpd2 




'■pp 
iVSK, 



pd.2 







2K. 



pp 



iVSKpd2 K^ J 



Ki -2iKzd, 
2iKzd, K^ 

K° 

/ Ki VSKpd, \ 

VSKpd, K^ 2Kpp 

Ki VSKpd, 

V 2Kpp VSKpd, K^ J 

I 



(37) 



(38) 



(39) 



(40) 



We can notice that Eq. 



has the same struc- ture as ([5^ . with two 2x2 degenerate sub-blocks hy- 
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bridized by a non-diagonal clement {Kpp in this case). 
This results in a splitting of the single-layer levels 
Epd2,+ {K) -^ Epd-2.+a{K),Epd2,+b{K), Epd2,~{K) -^ 
Epd2-a{K),Epd2-b{K). The two levels Epd2,+a{K), 
Epd2,+b{K), by looking at their orbital character, can 
be identified in DFT results in the small splitting of the 
(B) Epd2,+ {K) level, confirming once more the smallness 
of the interlayer Px/y-Px/y hopping. 

Less straightforward is the case of the 4x4 block 
Hpzdo2 i-^) where the hybridization term v2^pz niixes 
two different 2x2 sub-blocks, Hpdg and Hzd2- In this 
case, a mixing of the orbital character will result. Wc 
note, however, that the block Hpzdo2 i-^) appears twice 
in ([55)1 . so that each energy level will result double- 
degenerate, in particular the minimum (A) of the con- 
duction band at K. Note, however, that the negligible 
shift of such energy level in the DFT calculations with 
respect to the single-layer case is an indication that also 
the interlayer hopping element Kpz, between pz on one 
layer and Px, Py on the other one, is negligible. 



C. Q point 

An analytical insight on the electronic structure at the 
Q point was not available in single-layer systems and it 
would be thus even more complicate in the bulk case. A 
few important considerations, concerning the minimum 
(D), can be however drawn from the DFT results. In par- 
ticular, we note that in the single-layer case this energy 
level had a non- vanishing pz component. As we have 
seen above, the interlayer hopping between pz orbitals 
appears to be dominant with respect to the interlayer 
hopping between Px/y and Px/y and with respect to the 
mixed interlayer hopping Pz-Px/y We can thus expect a 
finite sizable splitting of the (D) level, containing a finite 
Pz component, with respect to the negligible energy shift 
of Epda.+ (A), which depends on the mixed interlayer 
process Kpz. 



V. MOMENTUM/ORBITAL SELECTIVE 

SPLITTING AND COMPARISON WITH DFT 

DATA 

In the previous section we have elucidated, using a 
tight-binding model, the orbital character of the band 
structure of M0S2 on the main high-symmetry points of 
the Brillouin zone. Wc have shown how a reliable min- 
imal model for the single-layer case needs to take into 
account at least the Px, Py orbitals of the S atoms in ad- 
dition to the Ad orbitals of Mo. A careful inspection of 
the electronic structure shows also that the band edges 
at the K point defining the direct band gap in the single- 
layer case are characterized not only by a chiral order of 
the d Mo orbitals, as experimentally observed, but also 
by an entangled chiral order of the minor component of 
the Px/y S orbitals. 



An important role is also played by the Pz orbitals 
of the S atoms. In single-layer systems, the Pz orbital 
character is particularly relevant in the (C) state, char- 
acterizing a secondary maximum in the valence band at 
the r point, and in the (D) state, which instead provides 
a secondary minimum in the conduction band at the Q 
point. 

The Pz component becomes crucial in multilayer com- 
pounds where a comparison with DFT results shows that 
the interlayer coupling is mainly driven by the Pz-Pz hop- 
ping whereas Px/y-Px/y^ Pz-Px/y are negligible. This re- 
sults in an orbital-selective and momentum-dependent 
interlayer splitting of the energy levels, being larger for 
the (C) and (D) states and neghgible for (A) and (B). 
This splitting is thus the fundamental mechanism respon- 
sible for the transition from a direct (A)-(B) gap in single- 
layer compounds to an indirect (C)-(D) gap in multi- 
layer systems. Controlling these processes is therefore of 
the highest importance for electronic applications. Note 
that such direct/indirect gap switch is discussed here in 
terms of the number of layers. On the other hand, the 
microscopical identification of such mechanism, which is 
essentially driven by the interlayer coupling, permits to 
understand on the physical ground the high sensitivity 
to pressure/strain effects, as well as to the temperature, 
via the lattice expansion. 

Finally, in order to show at a quantitative level how the 
orbital content determines the evolution of the electronic 
structure from single-layer to multilayer compounds, we 
have performed a fitting procedure to determine the 
tight-binding parameters that best reproduce the DFT 
bands within the model defined here. The task was 
divided in two steps: i) we first focus on the single- 
layer case to determine the relevant Slater-Koster intra- 
layer parameters in this case; ii) afterwards, keeping 
fixed the intralayer parameters, we determine the inter- 
layer parameters. To this purpose we employ a simplex 
mcthodiSSi to minimize a weighted mean square error 
/wMSE between the TB and DFT band energies, defined 
as 



/wMSE = 2_^ ^''' 



(k)[e™(k)-e?^T(k)]^ 



(41) 



where eP^'^(k) is the dispersion on the i-th band of 
the 11 band block under consideration, e™(k) the 
corresponding tight-binding description, and Wi(k) a 
band/momentum resolved weight which can be used to 
improve fitting over particular k-regions or over selected 
bands. In spite of many efforts, wc could not find a 
reliable fit for the whole electronic structure including 
the seven valence bands and the four lowest conduction 
bands ii£^ As our analysis and our main objective con- 
cerns the description of the valence and conduction bands 
that define the band gap of these systems, we focus on 
finding a set of parameters that describe properly these 
bands. Since both the lowest conduction and highest 
valence band belong to the electronic states with even 
z ^ —z symmetry, the fit was performed in the 6x6 
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FIG. 5: Comparison between the DFT band structure (black 
dots) and the best fit tight-binding model (red solid lines) for 
single-layer (top panel) and bulk M0S2 (bottom panel). 



orbital space defined by this symmetry. In addition, due 
to the degeneracy at the F point and to the band cross- 
ing along the F-M direction, the two conduction bands 
with even symmetry for z — >■ — z were considered in the 
fit. Additionally, wc give a larger weight to the (A)-(D) 
band edges in order to obtain a better description of the 
most important features of the band structure. 

Our best fit for the single-layer case is shown in the 
top panel of Fig. [5] (where only the TB bands with even 
symmetry z — > — z are shown), compared with the DFT 
bands, and the corresponding tight-binding parameters 
are listed in Table |lll Note that, due to the restriction 
of our fitting procedure to only some bands belonging to 
the block with even symmetry, the atomic crystal field 
Ai for the Mo orbitals d^z, dyz (not involved in the fit- 
ting procedure) results undetermined. The fit reported 
in Fig. [5] agrees in a qualitative way with the DFT re- 
sults, showing, in particular, a direct gap at the K point 
[(A) and (B) band edges] and secondary band edges for 
the valence and conduction bands lying at the F (C) and 
the Q point (D), respectively. 

Turning now to the bulk system, the further step of de- 
termining the interlayer hopping parameters Uppc, Upp-^, 
is facilitated by the strong indication, from the DFT anal- 
ysis, of a dominant role of the interlayer hopping between 
the Pz orbitals and a negligible role of the interlayer hop- 
ping between the Px/y orbitals. Focusing on the F point, 
these two different hopping processes are parametrized in 
terms of the corresponding interlayer parameters Tzz and 
F, 



and UppTT, by fixing the effective splitting of the Ezdo.+{^) 
level as in the DFT data. The values of Upp^ and t/pp^ 
found in this way are also reported in Table HIl and the 
resulting band structure in the lower panel of Fig. [SJ 
where only the TB bands with even symmetry z — >■ — z 
are shown. We stress that the intralayer hoppings are 
here taken from the fitting of the single-layer case. The 
agreement between the DFT and the tight-binding bands 
is also qualitatively good in this case. In particular, we 
would like to stress the momentum/orbital selective in- 
terlayer splitting of the bands, which is mainly concen- 
trated at the r point for the valence band and at the Q 
point for the conduction band. This yields to the crucial 
transition between a direct gap in single- layer M0S2, fo- 
cated at the K point, to an indirect gap F-Q in multilayer 
systems. 

On more quantitative grounds, we can see that, while 
the interlayer splitting of the condution level EzdQ,+{^) is 
easily reproduced, the corresponding splitting of the con- 
duction band at the Q point is somewhat underestimated 
in the tight-binding model (0.42 eV) as compared to the 
DFT data (1.36 cV). This slight discrepancy is probably 
due to the underestimation, in the tight-binding model, 
of thepz character of the conduction band at the Q point. 
As a matter of fact, the set of TB parameters reported 
in Table HIl gives at the Q point of the conduction band, 
for the single-layer case, only a 3.8% oi pz orbital char- 
acter, in comparison with the 11% found by the DFT 
calculations. It should be kept in mind, however, that 
the optimization of the tight-binding fitting parameters 
in such a large phase space (12 free parameters) is a quite 
complex and not univocal procedure, and other solutions 
are possible. In particular, a simple algebric analysis sug- 
gests that an alternative solution predicting 11% oi pz 



Crystal Fields 


Ao 
Ai 
A2 
Ap 
A. 


-1.016 

-2.529 
-0.780 
-7.740 


Intralayer Mo-S 


Vpda 

VpdTV 


-2.619 
-1.396 


Intralayer Mo-Mo 


Vddcr 
VddTT 
VddS 


-0.933 

-0.478 
-0.442 


Intralayer S-S 


Vppcr 

VppTT 


0.696 

0.278 


Interlayer S-S 


tVppTT 


-0.774 
0.123 



- pp, as discussed in Appendix |B] We can thus approxi- 
mate Tpp = 0, providing a constraint between Uppa and 
UppTt , and leaving thus only one effective independent fit- 
ting parameter: Tzz- We determine it, and hence Uppa- 



TABLE II: Tight-binding parameters for single-layer M0S2 
(Act, Va) as obtained by fitting the low energy conduction 
and valence bands. Also shown are the inter-layer hopping 
parameters Ua relevant for bulk M0S2. All hopping terms 
Va, Ua and crystal fields A^ are in units of eV. 
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character at the Q point would yield to a corresponding 
splitting of the order of 1.2 eV, in quantitative agreement 
with the DFT data. A more refined numerical search in 
the optimization of the tight-binding parameters, using 
global minimization techniques, might result in better 
comparison with the DFT results and further work along 
this line should be of great interest. 



We define thus: 



VI. CONCLUSIONS 

In this paper we have provided an analytic and reli- 
able description of the electronic properties of single-layer 
and multi-layer semiconduting transition-metal dichalco- 
genides in terms of a suitable tight-binding model. We 
have shown that the band structure of the multilayer 
compounds can be generated from the tight-binding 
model for the single-layer system by adding the few rele- 
vant interlayer hopping terms. The microscopic mech- 
anism for the transition between a direct-gap to an 
indirect-gap from single-layer to multi-layer compounds 
is thus explained in terms of a momentum/orbital selec- 
tive interlayer band splitting, where the orbital Pz com- 
ponent of the S atoms plays a central role. The present 
work provides with a suitable basis for the inclusion of 
many-body effects within the context of Quantum Field 
Theory and for the analysis of local strain effects related 
to the modulation of the Mo-S, Mo-Mo and S-S ligands. 
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Appendix A: Tight-binding Hamiltonian elements 

In this Appendix we provide an analyical expression, 
in terms of the Slater-Koster parameters, for the several 
intra-layer and inter-layer matrix elements that appear 
in th e Ha miltonian of the tight-binding model. Following 
Ref. Il04 it is convenient to introduce few quantities that 
account for the moment dispersion within the Brillouin 
zone, as functions of the reduced momentum variables 
^ = kxa/2, T] = y/3kya/2. 



Ci(C,?y) = 2 cos(C) cos(r//3) + cos(2?y/3) 

+i[2 cos(C) sin(?7/3) - sin(2r//3)], (Al) 
C2{^,V) = cos(^) cos(ry/3) - cos(277/3) 

+i [cos(0 sin(77/3) -t- sin(2ry/3)] , ( A2) 
C3(C,7y) = cos(Ocos(?7/3)-H2cos(27;/3) 

+i[cos(0 sm{i]/3) - 2 sin(27//3)], (A3) 

dii^,7j) = sin(77/3)-^cos(7,/3), (A4) 

/i(^,77) = cos(20 + 2cos(Ocos(?7), (A5) 

l2i^,v) = cos(20-cos(Ocos(7?), (A6) 

hi^^Tj) := 2cos(20+cos(e)cos(77). (A7) 
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1. Intra-layer hopping terms 



Following Ref. |104 the intralayer hopping terms Ha^p 
appearing in Eqs. (|4])-([8]) can be written as: 



H.. 



y/y 



H^ 



H. 



z^ /z^ 



H, 



X^ I X^ 



H. 



xy/xy 



H. 



xz/a 



H, 



yz/y-' 



H, 



/y 



H. 



^ z^ /xy 
^x'^ / xy 

TT 

^^xz/yz 



Hz2/x 



H. 



' !v 



H. 



H.. 



X^ / X 

x^/y 



H. 



xy I X 



H 



xy/y 
xy/z 



H. 



XZJ X 



H 



xz/y 



H. 



XZ I Z 
TT 

-^yz/x 






yz/y 

yz/z 









^ Ap + Eizhi^, V) + ^Eie cos(C) cos(?7), 
= Ap + Eiehi^, V) + 3-Ei5 cos(C) cos(?7), 
= A, + 2Ei6lii^,r,), 



A2 + E^M^, V) + 3^12 cos(0 cos(7?), 
A2 + Ei2h(.C, V) + 3^11 cos(0 cos(7?), 
Ai + EishiC, 77) + 3^14 cos(0 cos(77), 
Ai + EuhiC, 77) + 3^13 cos(0 cos(7?), 

-V3(£;i5-Si6)sin(e)sin(77), 

2Ewl2itv), 

-2^3^10 sin(0sin(7?), 

%/3(-Bii-Si2)sin(Osin(77), 

V3(^i4-£i3)sin(e)sin(?/), 

-2^/3i;lSin(Odl(e,r;), 
2E^C2{^,r,), 

-2V3(ii?5-i?3)sin(Odi(^,?7), 
-2i?3C3(C, ??) - 2*£;5 cos(Odi(C, r?), 
-2i?4C2(C,?7), 

-IezCsI^^v) - Q^EsCosiOdiitv), 

Hx^/x{i,v)i 
2V3Eisin{Odi{^,7j), 

^EeCsit v) + QiEj cos(Odi(^, v), 

2V3{^Ee~Er)smiOdiiC,v), 

-2V3Essm{^)di{^,Tj), 

Hxz/y{^,V), 

2ErC3{tv) + 2iEecos{0diiC,v), 
2EsC2i^,v), 



where 



El 



E2 = 



Vpda f sin^ - ^ cos^ j + VSVpdTT sin^ 1 



1 
2 

X cos 4>, 

-Vpda [ sin 



(A8) 



■ cos 



-VSVr 



pd-K COS 



X sm 0, 



-E3 = - 



Ea 



73 
2 

V3 



— ypda cos"* 4) + VpdTT COS 4) sin 



Vnda sin (j) cos'^ (j) — VpdTT sin (p cos^ 1 



^ pda ' 



E5 = - -VpdT, COS (f), 

3 

^6 = - -Vpd-^ sin (f), 



Ej 



E. 



1 
4 

X sin 
1 
2 
X cos 



V^Vpda COS^ <j> ~ Vpdn{l - 2 cos^ , 

0, 

V^Vpda sin^ (j) - VpdTT {I - 2 sin^ q 



1 3 

E() = -Vdda + -rVddS, 

V3 



(A9) 
(AlO) 

(Ml) 

(A12) 
(A13) 

(A14) 

(A15) 

(A16) 

(A17) 

(A18) 

(A19) 
(A20) 
(A21) 

(A22) 
(A23) 

Here the angle (j) characterize the structure of the unit 
cell of the compound and it is determined by purely ge- 
ometric reasons (see Fig. [T]).. For the ideal trigonal 
prism structure, neglecting the marginal deviations from 
it in real systems, we have (p = arccos[-\/4/7], so that 
cos (j) ~ \/4:/7 and sin (j) ~ \/oj7. 

With these expressions, taking into account also the 
further changes of basis, the Hamiltonian at the F point 
can be divided in sub-blocks as: 

Hzdoir) 

i?E(r) = I HpdAT) I ,(A24) 

Hpd,{r) 



Elo 


— . [Vdda Vdds] , 


Ell 


3 1 

— -jVdda + -rVddS, 

4 4 


E12 


= VddTT, 


Eu 


= VddT,, 


El4 


= Vdds, 


Ei5 


= Vppa , 


Eie 


= VppTT ■ 



Hoir) = 



HpdA^) 

HpdAr) 

OF, 



(A25) 
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where 



H.doiT) 



To -J^Tzdo 



2 do 



\/2r 



ffp..(r) - ( ^ p; 



Pd2 



P 



(A26) 



(A27) 



ffp<i.(Al = 



HpdAK) 



K2 iVSK, 



pd2 



, (A45) 



(^aU1§''|' <*^^' 



Hpdi (r) 






(A28) 



The parameters Fq can be viewed as "molecular" en- 
ergy levels, and the quantities Ta^p as hybridization pa- 
rameters. Their explicit expressions read: 



= Ai + 3[Ei3 + Eu], 

= A2+3[^11+Si2], 



r^ = r,. 



Vr 



PPTTJ 



pO 

p 

pE 
T-O 



r — v^ 









ppcrj 

= Ap + 3[Eir_, + Ei6], 



H- 



S.z'^-r'^/z 



(r) - iE2, 



i pd2 



-f^2;2-y2/y(^) 



iJ, 



xy /x 



(r) 



= -2[3S3 + £;5], 

= -f^2;2/2;(r) = Hyz/y{T) 

- 2[3^7 + ^6]. 



(A29) 

(A30) 

(A31) 
(A32) 

(A33) 
(A34) 
(A35) 

(A36) 
(A37) 
(A38) 

(A39) 

(A40) 



At the K point, in the proper basis described in the 
main text, we can write the even and odd blocks of the 
Hamiltonian as: 



H^{K) 



He 




where 



^^'^o(^) = I 2.S„ 



Hzd2{K) 



K2 
2Kzd2 



'2iKpdo 



2K,a2 



(A41) 



(A42) 



(A43) 



(A44) 



HzdAK) 



2iK,r. 



-2iK,d^ 
K? 



The parameters Ka, Ka^p read here: 



Ko 
Kl 

Ko 



p 
K^ 

Kl: 

K? 
K„ 



H,2/,2{K) = An-3E<,, 

HXZ/XZ{K) = Hy^/y^{K) 

Ai--[Ei3 + Ei4], 



Hxy/xy{K) 

2^ 
-K, 



A2 
K 



,[Eii 



- Hx^/x'^iK) 
\- E12], 



p ^ ^'ppTT ; 
^^p y ppiT ^ 

K, - K, 



= K, 
= H. 



*ppij 1 



K, = 

Kpda = 

K,,d2 = 

Kpd.2 = 

K„d, = 



./X{K) = Hy/yiK) 

3 

^p ^ ^[^15 + Eie], 

H,/,{K) = A,~3Eie, 

Hsz^-r-^/yiK) = iH3^2_r2/x{K) 

-3^1, 

Hx2-y-^/z{K) = iHxy/z{K) 

3^4, 

H 



x2^y2/y{K) = -Hxy/xiK) 

-iHx2_y2/^{K) = -iHxy/y{K) 

[E5 - 3E3] , 

HXZ/X{K) = -Hy^/y{K) 

= iHxz/y{K) = iHyz/x{K) 
= [Ee^3E,], 
Kzd, = Hy^/ziK) = iHx,/,{K) 
= -SEs. 



^pd 



(A47) 



(A48) 
(A49) 

(A50) 

(A51) 
(A52) 
(A53) 
(A54) 

(A55) 
(A56) 

(A57) 

(ASS) 

(A59) 

(A60) 
(A61) 



2. Inter-layer hopping terms 

Inter-layer hopping is ruled by the Slatcr-Koster pa- 
rameters t/ppcr, UppTT describing hopping between S-3p 
orbitals belonging to different layers. 

In terms of the reduced momentum variables ^ = 



kxa/2, 7] = \/3kya/2, we have thus: 
I./A^,V) = ^[EigC3i^,-fj) +t3E,r cos ^dii^, ^4^^) 

ly/y {^,V) = l[El7C3i^,-v)+ ^^E.g COS ^dl {^ ^^P?) 



i./.M,v) 


= EisC,{^,^v), 


(A64) 


Ix/y{£..V) 


= ^[i?i7-i?i9]sin^di(e,-?y), 


(A65) 


h/zii.r,) 


= -\/3£'2osm^(ii(^, 77), 


(A66) 


ly/zi^iV) 


= -E20C2{^,~V), 


(A67) 


Iz/z{^,V) 


= ii;i8Ci(e,-r,), 


(A68) 


where 








£'17 = Uppa cos^ /^ + UppTT sin^ /3, 


(A69) 




£^18 = C/ppff sin^ /3 + C/pp^ cos^ /3, 


(A70) 




£19 = C^ppTT, 


(A71) 




£20 = [C/pp<T - C/pp^] cos /3 sin /3. 


(A72) 



Here /3 is the angle between the hne connecting the 
two S atoms with respect to the S planes (see Fig. [T]). 
Denoting w the distance between the two S-planes, we 
have: 
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Appendix B: Decomposition of the Hamiltonian in 
sub-blocks at high-symmetry points 



In this Appenddix we summarize the different unitary 
transformations that permit to decomposed at special 
high-symmetry points the higher rank Hamiltonian ma- 
trix in smaller sub-locks. In all the cases we treat in a 
separate way the "even" and "odd" blocks, namely elec- 
tronix states with even and odd symmetry with respect 
to the z — >■ — z inversion. 



cos/3 = 
sin/3 = 



ya? + 3uP 



V^v 



(A73) 
(A74) 



Using typical values for bulk M0S2, a = 3.16 A, and 
w = 2.975 A, we get cos/3 = 0.523 and sin/? = 0.852. 
At the high-symmetry points F, K, we have thus: 



Fpp 


= 4/.(F)=/,/,(F) 






3 


(A75) 


Tzz 


- Iz/z{T) 






= 3£'i8, 


(A76) 


Kpp 


= 4/,(/v) = -ly/y{K) 

= -iIx/y{K) = -iIy/^{K) 






= 7 [£"19 ^ £17] 1 
4 


(A77) 



Kp. = Iy,AK)^h,y{K) 

= -iI^/AK) = -tlz/AK) 
3 



-E- 



20- 



(A78) 



1. Single-layer 



a. F point 



In the Hilbert space defined by the vector basis (^J, in 
Eq. ^ , the even and odd blocks of the Hamiltonian can 
be written respectively as: 



17 



^E(r) 



( To 














y2r,d„ \ 





Tz 








V^pd^ 











r2 


v2rpd2 














VZFpd^ 


pE 

V 











\/2rpd2 








pE 
V 





V \/2r,do 














r? / 



(Bl) 



and 



i?o(r) = 



Ti 





vsrpdi 





\ 





Ti 





VZFprfj 





vsrpdi 





pO 











vsrpdi 





pO 

















r^/ 



(B2) 



The division in sub-blocks is already evident in Eqs. 
Bip - (|B2p . They can be further ordered using the basis 



h. K point 



4>l - 

where 



'^k.zda^'Pk' 



,iy' Vk,pd2,x^ 'Pk.pdi,x^Vk,pdi.y^ '?'fe,y?W 



k.zdn 



k,pd2,y 

P 
k,pd2.,x 

At _ Mt ^t 



\"'k,x'2-y^^Pk,y,s)-: 



\"'k,xy^Pk,x,s)^ 
^k,pdi,x ~ \"'k,xz^Pk,x,A)^ 
°k,pdi,y ~ \"'k,yz^Pk,y,A)^ 



(B4) 
(B5) 
(B6) 
(B7) 
(B8) 



<^L = K.z.s)- 


(B9) 


In this basis we eet Eos. ()ll|)-p^. where 






(BIO) 


n , (v\ ~ ( ^2 V2rpd2 \ 


(Bll) 


A- /px _ / Ti \/2rprfi \ 


(B12) 

1 



In the basis defined by the Hilbert vector (t>i, the 



^fc' 



Hamiltonian at the K point reads, for the even and odd 
blocks, respectively: 



Mk) 



Ko 








-iV2Kpd„ 


V^Kpdo 








K2 





iV2Kpd2 


V2Kpd2 


V2Kzd2 








K2 


~V2Kpd2 


iV^Kpd2 


-iV2K,d 


iV^Kpdo 


-iV2Kpd2 


-V2Kpd2 


p 








V2Kpdo 


V^Kpd2 


~iV2Kpd2 





K^ 








V2K,d2 


iV2K,d2 








Kf 



(B13) 
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HoiK) = 



Ki ~iV2Kpd, -V2Kpd, V2Kzd, 

V2Kpd, iV2Kpd, K° 

iV2Kpd, -V2Kpd, K^ 

\iV2K,d, V2K,d, K^ / 



(B14) 



r 



In order to decoplcd the Hamiltoniam, it is convenient 
to introduce the chiral basis defined by the vector V'l in 
([T7]). In this Hilbert space we have thus: 

HpdAK) 

He{K) = \ H^dAK) I (B15) 

Hpd,{K) 

Hpd, [K) 

Ho = \ H.dAK) 0_ I , (B16) 



KO 



where 



■0 ^^p 



HpdM) - y^2Kpd, le 
fj , (K\ - ( -^2 '^Kpd; 



(B17) 
(B18) 



K2 iVSKpd, 

^pd2 ^p 



Hv^AK) = ( ,.;^t, . --T' ), (B19) 



iV^Kpd, Kl 



fr (K\ - ( ^'^ -i^Kpd, 



(B20) 
(B21) 



2. Bulk system 

The general structure of the tight-binding Hamihonian 
Hhuik for the bulk system, using the basis defined in p3| , 
is provided in Eqs. (P5)) - (PT|) . where we also remind the 
symmetry property (|25p that related the matrix elements 
of iJ2 to Hi. 

As mentioned in the main text, for fc^ = the band 
structure can be still divided in two independent blocks 
with even and odd symmetry with respect to the trans- 



formation z — > 
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Further simplicication are encoun- 



tered at the high-symmetry points T and K. 



a. r point 

We first notice that at the F point the relation ([^5)) 
does not play any role, i.e. ^2(F) = -ffi(F), where ^i(F) 
is defined by Eqs. (|T0| - (fT2|) in the main text. 

The Hamiltonian is thus completely determined by the 



interlayer hopping matrix / that at the F point reads: 



/(F) 



pp 











i^pp 











l\z 



(B22) 



A convenient basis to decoupled the Hamiltonian in 
smaller subblocks is thus: 



where 



\f 



k V fe,2do' k,pd2,y^ k,pd2.x^ k,pdi,x^ k,pdi,y^ ^iPi^^ 



st 



st 



$ 



t 

k.zdn 



(4.3z^-rM'PL,A,l'4.3.2-r2 2'PL,i?^4) 



k,pd2 ,x 

t 



J 



J 



$t 



k,pd2,y ~ y"-k,x^-y^,l^Pk,y,S,l^"-k,x^~y^,2^Pk,y,S,2\l'^^'^) 

(B26) 
(B27) 
(B28) 
(B29) 



\"'k,xyS^Pk,x,S,l^"'k,xy,2^Pk,x,S,2)' 
J ^t ^t 



k,pdi,x ~ \"'k,xz,l'Pk,x,A,l^"'k,xz,2^Pk,x,A,2)' 
k,pdi,y ~ \"'k,yz,l'Pk,y,A,l'"-k,yz,2TPk,y,A,2)' 
^ yPk,z,S,l'Pk,z,S,l'^' 



H.. 



The resulting total Hamiltonian can be written as: 

HE,buik(F) ^ _ V (B30) 

iio,buik(F) y ' ^ ' 



-ffbulk(F) 



where 

^E,bulk(F) = 
^0,bulk(F) = 

and where 



-ffzdo,bulk(F) 

-ffpd,,buik(F) ^ (B^l,) 

-ffpd2.bulk(F) 

-ffpdi,bulk(F) 

-ffpdi,buik(F) ($32) 

i?.,buik(F) 



H 



2:tio,bulk 



/ Fo \/2Tzdo 

V2r,do rf 







H, 



pd2 jbulk 



H, 



pd\ ,bulk 



H, 



bulk = 



V 

v2Fpd2 


V 

v2Fpdi 



Tzz F? 



V2F, 



pE 
P 





pd2 






To 

vZFzrfo 




F2 

^ pp V ^-L pd2 

v2Fpdi 
FO 

Fi 

r 





^ z 





pE 
P 








v^f:./^^^) 



V2Vpd, fo ; 

(B36) 



b. K point 

The treatment of the bulk Hamihonian at the K point, 
in order to get a matrix elearly divided in bloeks, is a bit 
less straighforward than at the F point. 

We first notice that the interlayer matrix, in the basis 
!>{., reads: 



I[K) 



K 



pp 



iK, 



pp 



^^pp ^pp 



iK, 



pz 



Kp, 




(B37) 



* 
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We then ridefine the orbitals d], 

'"■k,yz,2' Pk,y,a,2 



d^ 



k,yz,2 k,yz,2 

Pk,v,a,2 = -P'k,v,a,2 (" =A,S), in 

order to get, according with (US]), ^2(r) = Hi{r). 



J 



nt 



Following what done for the single layer, we can also 
introduce here a chiral basis. After a further rearrange- 
ment of the vector elements, we define thus the conve- 
nient Hilbert space as: 



\^ k,pzdo2,L' ^ k,pzdo2,R' ^ k,pd2.E^ ^ k,pzdi,R 



' ^k,pzdi,L' ^k,pdi,0)^ 



(B38) 



where 



k,pzdo2,L 



* 



t 

k.pzdo2,B, 



* 



* 



t 

k,pzdi .R 



* 



k,pzdi ,R 
^k,pdi,0 



('^fc,3z2-r2,l'Pfc,L,S,l''^fc,i?,2'Pfc,z,A,2('v39) 
^"■k.az^ -r^ ,2^ Pk,B.,S,2^ "'k,L,l^ Pk,z,A,lW^^) 

{4ma^pIz,s,i^pIr,a,2)^ (B42) 

(4,l,2'PL,s,2'PLl,A4)' (B43) 

(4.La'PLfl,Aa'4,i?,2'PU,A,2J- (B44) 



In this basis, the Hamiltonian can be once more written 



as: 



Hhuik{K) 



H^{K) 
Ho{K) 



(B45) 



where H-e{K), Ho{K) are defined in Eqs. dMI-dini) of 
the main text. 
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